function res = central_der1_onestep(u,N,x_max)
h = x_max/N;
u = fill_ghostcell_period(u,N,1);
res = (u(3:N+3)-u(1:N+1))/(2*h); 